Discrete quantum gravity: continuum limit and the problem of state doubling 
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It is shown that in the theory of discrete quantum gravity defined on the irregular "breathing" 
lattice, if the macroscopic continuum phase is realized, the phenomenon of state doubling (even if 
it exists formally at kinematic level) actually is absent at experimentally accessible energies. 
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1. In the recent my works one more version of the dis- 
crete quantum gravity has been introduced into consid- 
eration 0-0. This version of discrete quantum gravity 
manifests "naive" macroscopic continuum limit. But it 
does not mean that in the theory the macroscopic con- 
tinuum phase indeed is realized dynamically. Moreover, 
some phenomena and problems are well known which pre- 
vent microscopic discrete quantum gravity from develop- 
ing into macroscopic continuum phase with the proper- 
ties close to the real one. For example, the phenomenon 
of state doubling (Wilson doubling) and the cosmological 
constant problem. 

In the present and the subsequent letters I study the 
state doubling phenomenon and the cosmological con- 
stant problem under the assumption that the macro- 
scopic continuum phase is realized dynamically in the 
theory of discrete gravity of the type |£g. I think that 
this is a very strong assumption: for realization of this 
assumption the microscopic theory must be determined 
very specially and, perhaps, uniquely (the correct set 
of fundamental fields and their interaction, and so on). 
Though all that is unknown, I apply a sort of phenomeno- 
logical approach and study the mentioned problems in 
the macroscopic phase. It seems to me that the knowl- 
edge of some qualitative properties of the considered 
theories in the macroscopic phase can promote general 
progress of this branch of quantum field theory. 

In the present letter I study the possible consequences 
due to the phenomenon of state doubling. The problem 
of graviton states doubling in discrete gravity was con- 
sidered earlier in ■ The fermion state doubling (proper 
Wilson doubling, see, for example 0) was investigated 
in 0-0. For simplicity I study here the Dirac field on 
the lattice. But it is clear from the consideration given 
below that the conclusion remains valid with respect to 
all other fields. 

At first sight it seems that the theories in which the 
state doubling phenomenon is present (I call them as 
anomalous theories) must be very different from the usual 
theories (without state doubling) in the continuum phase. 



Indeed, in anomalous theories the transitions from the 
states with normal modes into the states with anomalous 
modes (and vice versa) exist. But the such transitions are 
not registered at experimentally accessible energies. 

The main result of the paper is the behaviour in the 
x-space of the propagator describing the propagation of 
anomalous modes 
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Here and below a is a certain constant that has the length 
dimension and is of the order of the effective lattice step 
and A is the number of the order of unit. But the propa- 
gators describing the propagation of normal modes have 
the conformal behaviour at the enough large energies. 
Thus, the propagators of anomalous modes decrease too 
fast as compared with the propagators of normal modes. 
Therefore the transitions between the states of normal 
and anomalous modes are strongly depressed at the scales 
much greater than a. The situation is similar to the fol- 
lowing one: the electron loop corrections to the external 
smooth electromagnetic field die out at the scales much 
greater than 1/m (m is the electron mass) . 

2. It is necessary to write out some formulae from the 
work 0. Here the notations are completely identical to 
that in 0. 

Let ^ be a 4-dimensional simplicial complex such that 
the 3-dimensional complex 6 = d & has the topology of 
3-sphere S 3 . To each vertex a, S ^, the Dirac spinors ipi 
and ip i belonging to the complex Grassman algebra are 
assigned. To each oriented edge cuaj G .ft, an element of 
the group Spin(A) 
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[7 a ,7 b ], (2) 



and also an element = 7°, such that 
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are assigned. The notations ipAi> ipAi, &Aij> ^Aij and 
so on indicate that the edge Xf- — cuaj belongs to 4- 
simplex with index A. Let aAi, aAj, aAk, clai, and 
a. Am be all five vertices of a 4-simplex with index A and 
£Aijkim = ±1 depending on whether the order of vertices 
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o-Ai ciAj o-Ak a>Al a Am defines the positive or negative ori- 
entation of this 4-simplex. The Euclidean action of the 
theory has the form 

I = 5x"24 Hi,j,k,l,m £ Ai]klm tr 7 5 X (4) 

X ^ 21'^ ^Ami^Aij ^Ajm^-Amh^Aml~\~ 
+ igJ a (4>Ai'J a ^Aij1pAj ~ ^ A ftAjil a i>Ai) eArnjeAmkeAml} 

The oriented volume of a 4-simplex with vertexes 
a-Ai, a A], a Ak , a Al , and a Am is equal to 
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The quantity Ifj — J2 a =i( e ij) 2 1S interpreted as the 
square of the length of the edge ajOj. 

The partition function Z for a discrete Euclidean grav- 
ity is defined as follows: 
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Here, the vertices and edges are enumerated by indices 
V and £ , and the corresponding variables are denoted by 
"0v i ^£ ; respectively, d is the Haar measure on the 
group Spin(4), and 

de £ = ]Jde£, di/v dVv = J| dipvu dipw ■ (?) 



The index v enumerates individual components of the 
spinors ipy and 4>v 

Let's denote by X a 4-dimensional smooth manifold 
with topology of the complex ^. Consider a set of 
maps {g} from geometrical realization of the complex 
^ onto manifold X which are not necessarily one-one 
maps. For a given local coordinates x^, \l = 1, 2, 3, 4 
a map g defines the coordinates of images of vertexes 
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g^iflAi)- Define the four differentials 

j, i = 1, ...,4. 
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X g(Ai) ~ X g(Aj) ' 
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Suppose the smooth fields w'^i), e^(x), ip(x), ip(x) are 
defined on the manifold X. Then we can define the dis- 
crete lattice variables according to the relations 
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^(lj(Am))dlJ (Aml) =e^ mI , lp(x g (Ai)) = *pA t - (9) 

On the contrary, the discrete lattice variables in the 
right hand sides of Eqs. © define the values of the 
fields on the images of vertexes of the complex. It 



is clear [j| that for the discrete lattice variables which 
change enough smoothly along the complex we obtain 
the enough smooth fields. Moreover, in this case the dis- 
crete action transforms to the well known continuum 
action 

I = Je abcd {-jrR ab /\e c Ae d + 
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uj ab = LU^fdX^, 
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I emphasize that we obtain the action 11(J|I only if the 
lowest derivatives of the fields are taken into account. 
It is important that in this case the information on the 
structure of the complex is lost. This is incorrectly if 
the highest derivatives of the fields are also taken into 
account. In the work |5( the arguments have been given 
that the quasi-classical phase at the same time is the 
macroscopic continuum phase with long correlations and 
hence, also, the phase in which the highest derivatives 
of the fields can be ignored. In this phase the partition 
function is saturated by normal (smooth) modes but not 
by anomal modes responsible for Wilson state doubling. 

It is important that in the quasi-classical phase the 
universe wave function does not depend on the discrete 
variables e Ai j in a wide diapason. This statement is true 
in the same context as the highest derivatives of the fields 
can be ignored. Indeed, in the quasi-classical phase the 
action H10J) as well as the universe wave function depend 
on the fields u>^ b (x), e^(x),ip(x),t/j(x) which are present 
in the left hand side of Eqs. ©■ Thus, fixing these 
fields and varying the maps g or, cquivalently, the images 
of vertexes and, hence, the differentials dx^ A j^, 

one can vary the discrete variables e A mi in the right hand 
side of Eqs. © in a wide region. 

To clarify the situation, let's consider the case when 
geometry of the space-time is flat or almost flat. In the 
flat case one can take 
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Here, the sum in the parentheses is taken over any closed 
path formed by 1-simpliccs. Equations (|llf> indicates 
that the curvature and torsion are equal to zero. Thus, 
geometrical realization of the complex & is in the four- 
dimensional Euclidean space, e?- being the components of 
the vector in a certain orthogonal basis in this space, and 
if Rf is the radius- vector of vertex a i: then efj = Rj—Rf. 
In this case one can take e°(a;) = <5°. It is evident that 
Eqs. are the only restrictions for the variables efj. 

3. Let us write out the equation for the eigenmodes 
of the discrete Dirac operator in partial case This 
operator is obtained by varying action Q with respect 
to the variable ip i . In the four-dimensional case, we have 
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(see 0) 



where 
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Here, the index A(i,j) enumerates all the 4-simplexes 
containing the edge aj<2j, the index enumerates all 
the vertices o,jU) neighboring the vertex a; @, and w,; is 
the sum of oriented 4- volumes of all 4-simplexes contain- 
ing the vertex a^. There exists the relations 

XX 4 =4^*, E^= ' ( 14 ) 

that provide the transformation of discrete equation 1 (12(1 
in the long-wavelength limit to the continuum equation 



i^ a d a tp = eijj. 



Indeed, this is easily verified by taking into account the 
relation 



(16) 



which is almost true in the long-wavelength limit, i.e., for 
conventional modes. Here, x a are the Cartesian coordi- 
nates in the same orthonormal basis in which the compo- 
nents of the vectors efj are specified. It is again seen that 
information on the positions of vertices of the complex 
is completely lost in the long-wavelength limit; i.e., con- 
ventional long-wavelength modes lose information on the 
lattice structure. The eigenvalues of the modes of Dirac 
operator (|15fl are e = ±|fc|, where k a is the wavevector 
of a mode. Thus, the eigenvalues of the usual modes can 
be arbitrarily small. 

The main result of works Q, Q is that the theory 
under consideration involves two types of lattices. On the 
lattices of one type, fermion doubling occurs, whereas it 
is absent on the lattices of the other type. 

According to the definition, the lattice allows fermion 
state doubling if the discrete Dirac operator given by 
Eq. Ijl 2(1 allows two types of modes with eigenvalues ap- 
proaching zerojldl| . The qualitative difference between 
them is as follows. The normalized modes of the first 
type, or normal modes, satisfy the conditions 



I AT 



-1/2 



where di and aj are the neigboring vertices, e is the mode 
eigenvalue, and N is the number of the vertices of the 
complex. The normalized modes of the second type, or 



anomalous modes, for certain neighboring pairs of ver- 
tices whose number is comparable with the total vertices 
number of the complex satisfy the conditions 
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Relations 1(170 mean that anomalous modes generally 
change in spurts from a vertex to a neighboring vertex. 
Hence, the derivatives d a ip A of anomalous modes also 
generally change in spurts from a vertex to a neighbor- 
ing vertex. Otherwise, Eqs. 1(12(1 and 1(14(1 would provide 
the equation ij a d a i/j A = e^" 4 , where the left-hand side 
is continuous, while the right-hand side is discontinuous. 
Moreover, the derivatives of anomalous modes in various 
directions are incommensurable. 

Further, we will show that the effective equation in the 
continuum limit that describes anomalous modes has the 
form 



ia a (x)d a g A {x) = eg A {x), 



(18) 



where a a {x) are random functions in the sense indicated 
(15) below. Indeed, let t/O , s — 1, 2, ... be a complete set 
of zeroth anomalous modes [llj. Any linear combination 
of zeroth modes is evidently also a zeroth mode. For this 
reason, we seek a soft anomalous mode "growing" from 
zeroth anomalous modes in the form tp A — g A jip^j°\ 
where the numerical field g A j is slowly varying. To derive 
an equation for the field g A j , we substitute the expression 
for a soft anomalous mode into Eq. 1 (12(1 : 
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The first equality in Eq. 1 (19(1 follows from the fact that 
4>fj°\s = 1, 2, ... are zeroth modes. According to the 
second equality in Eq. I(19|). the eigenvalue e for slowly 
varying fields g^ can be arbitrarily small. In this case, 
Eq. 1(19(1 is reduced to the form 
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The multiplication of this equation by the matrix 
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^• A(0 V^ (0) 



, reduces it to Eq. |JT%)|. where 
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Thus, the operator on the right-hand side of Eq. i|18|) 
is a first-order differential operator with variable matrix 
coefficients that acts on a multicomponent (or vector) 
function g s (x), s = 1, 2, . . .. 

Quantities H2U|) are irregular functions that depend 
strongly and locally on the positions of the vertices of 
the complex. This fact is demonstrated by the evident 
formulas in the two-dimensional theory (see Q). I give 
here only qualitative description of the quantities l|20(l 
and their properties. 

As it is known, the zeroth anomalous modes look like 
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where n s and n s (j) are whole numbers, n s > 1, n s (j) = 
0, ±1, . . . , ±(n s - 1), and n s (ji) ^ n s (j 2 ) for the neigh- 
boring vertices aj 1 and a j2 . From here it follows that 



, A(0) , ,A(0) 
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(22) 



for the neighboring vertices aj 1 and a j2 . In consequence 
of Eqs. (|21|l and 12211 the sums in the second square 
brackets in (|20|l become very irregular [l^. Even the 
signs of these sums and the oriented volume Vi in 114|) 
are not correlated. 

Thus, the components a a (x) in Eq. l|18f) are discontin- 
uous functions of x and they depend strongly and locally 
on the integration variables {eg} in the integral ©. It 
is important that the variables {e^} in integral 10 in 
the quasi-classical phase can independently vary over a 
wide range without a change in the action (see the text 
between Eqs. {THl and {TTJ). 

Further I assume that 



a a {x) =[3 a {x) + \p a , 



(23) 



where p® s , is the constant matrix, A is the small numer- 
ical parameter, and all odd powers of f3(x) in the inte- 
grand in Eq. © are equal to zero: 

(\0( Xl ) ...P(x 2n +i)\) e£ =0, 71 = 0,1,.... (24) 

Here the matrix product of matrices (3 is averaged. The 
points x\, ... , X2n+i can partially or completely coin- 
cide with each other. Averaging in Eq. I|24|) is performed 
by means of functional integral 10 • The subscript {eg} 
indicates that integration is performed only with respect 



to the variables {eg} 0. We also assume that the pa- 
rameter A is small so that the expansion in this parameter 
is meaningful. The last assumption is in agreement with 
said above. 

When the problem of the S matrix (which is meaning- 
ful only in the continuum limit) is solved in the theory 
under consideration, all S matrix elements mast be av- 
eraged over the variables {eg} before the calculation of 
probabilities and cross sections. Therefore, if the vertices 
describing the interaction are universal (i.e., independent 
of the microstructure of the lattice), the propagators of 
the matter fields in such a theory must be averaged over 
the variables {eg}. Indeed, in this case, the diagram- 
matic technique can be obtained in this case as a result 
of the expansion in a functional differential operator act- 
ing on the transition amolitudc for the matter fields in the 
quadratic (" free" ) approximation against the background 
of the external field sources. In this case, the square tran- 
sition amplitude must be averaged with a weight over the 
variables {eg } before the expansion in interaction. 

So, it is necessary to obtain the propagator 



■Mo) 



(25) 



that describes the propagation of anomalous modes. 
Thus, the problem is reduced to the study of the Green's 
function that corresponds to the operator on the left- 
hand side of Eq. fl%| and is averaged over the variables 
{eg} by taking into account conditions l|24(l : Il3 |: 



<(-^a a )- 1 ) e 



(26) 



4. We first consider the case A = 0. To perform the 
necessary averaging in Eq. I|25(l . we use the well-known 
formula 
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This formula implies that t / or r _1 ^ ±oo. Using 
Eq. H27|) . w e represent operator l|26ll (for the case A = 
before averaging over the variables {eg}) in the form 



(x\(-i(3 a d a ) X \y) 



1 



As (a 



e sf3 a d a _ e -sf3 a d a 



\y)- (28) 



This relation is meaningful for x =/= y, because its left- 
hand side is finite in this case. However, representation 
(|28|l is meaningless for x = y. Then, let us average the 
terms in the square brackets in Eq. I|28l) over the variables 
{eg}. In this case, it is only important that, owing to Eq. 
(|24|l . the result of such averaging of each of these terms is 
a function of (±s) 2 = s 2 but not of (±s). Therefore, the 
terms in the square brackets in Eq. I|28|l cancel each other 
after the averaging over the variables {eg} and thereby 
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— > — + > x > + ••• 
FIG. 1: 

the propagator describing the propagation of anomalous 
modes is proportional to the <5 function of the variable x: 

Ux\{-i^d^- X \y)\ ~a8W(x-v)- ( 2 9) 

I remind the reader that a is a certain constant that has 
the length dimension and is of the order of the effective 
lattice step. 

Now let's estimate quantity (|26|) for nonzero A by tak- 
ing into account that it can be expanded in the parame- 
ter A. This problem is easily solved by using field theory 
methods. Figure 1 shows the sum of the diagrams that 
corresponds to quantity 12tj|) . The solid line and cross 
correspond to unperturbed propagator (|29[l and the op- 
erator i\p a d a , respectively. Then using the conventional 
diagrammatic technique rules, one can represent quan- 
tity H'26l) in the form of the series of diagrams shown in 
Fig. 1. 

As a result, the S function of the free propagator is 
"smeared" and the desired propagator describing the 
propagation of anomalous modes appears to be exponen- 
tially decreasing in the x space according to the Eq. . 

At the same time, according to Eq. I|15(l . the free prop- 
agators describing the propagation of normal modes be- 
have according to the power law in x and momentum 
space. 

5. We briefly summarize the conclusions as follows. 

(i) On complexes that formally allow state doubling 
(of any fields), this phenomenon is really absent because 
normal and anomalous modes propagate differently and 
separately in spacetime. Normal modes (in the long- 
wavelength limit) have definite energy and momentum. 
On the contrary, anomalous modes cannot have definite 
energy and momentum and their propagators decrease 
exponentially in spacetime at scales comparable with the 
characteristic lattice scale. 

(ii) Let's give some general considerations concerning 
fcrmion (Wilson) doubling due to the above result. 

Let us assume that one Weyl field is introduced on the 
lattice by introducing the projection operator (1/2)(1 ± 
7 5 ) into the fermion action in Eq. |@}. Since the lat- 
tice fermion action is invariant under global 7 s trans- 
formations and the lattice fermion measure is invariant 
under local 7 s transformations, the total fermion cur- 
rent is strictly conserved. However, this does not mean 
that each of the currents of normal and anomalous modes 
is conserved separately. On the contrary, it has been 



well known for a long time that, when S matrix ele- 
ments are calculated by using causal Feynman propa- 
gators, the Weyl-field current is not conserved due to the 
gauge anomaly. This phenomenon is interpreted in the 
theory under consideration as the mutual transfer of the 
axial charge between normal and anomalous modes. Here 
it is shown that this phenomenon can occur only on the 
lattice scales. 

I also note that the theory under consideration does 
not exclude the case where the axial chargees of normal 
and anomalous modes are conserved separately (see 3]). 
This means that the axial anomaly is absent. Such a 
regime can be realized only in problems where the use 
of Feynman propagators is incorrect, for example, at the 
universe inflation stage when the problem of the S matrix 
is meaningless. 
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